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1. Introduction 

We report on our analytic studies of the renormalization properties of Borici-Creutz [JXJ, ||, 
[3|, and Karsten-Wilczek [BL ||] fermions (see [[7L |j| |9|], and references therein), two particular 
realizations of minimally doubled fermions. 1 These actions preserve an exact chiral symmetry 
for a degenerate doublet of quarks, and at the same time they remain strictly local, so that they are 
much cheaper to simulate than Ginsparg-Wilson fermions. They could then become a cost-effective 
realization of chiral symmetry at nonzero lattice spacing. This £/(l)<g)£/(l) chiral symmetry, which 
is of the same form as in the continuum, protects the quark mass from additive renormalization. As 
we have also verified at one loop, the renormalization of the quark mass has the same form as, say, 
overlap or staggered fermions. 

It is noteworthy that using minimally doubled fermions one can construct a conserved axial 
current which has a simple expression, involving only nearest-neighbour sites (see Section |6|). 
These actions are then among the very few lattice discretizations which provide a simple (ultralocal) 
expression for a conserved axial current. 

It is natural to compare these realizations of minimally doubled fermions with staggered 
fermions, which preserve the same U{\)®U{\) chiral symmetry and are also ultralocal and com- 
parably cheap. The advantage of Borici-Creutz and Karsten-Wilczek fermions is that they contain 
2 flavours instead of 4, and thus they do not require any uncontrolled extrapolation to 2 physical 



light flavours [Jllj, |12J. Moreover, the construction of fermionic operators is much easier than for 
staggered fermions, where there is also a complicated intertwining of spin and flavour. Minimally 
doubled actions look then ideal for Nf = 2 simulations. 2 



2. Actions 

The free Dirac operator of Borici-Creutz fermions is given in momentum space by 

d (p) = *£(y M sin/? M + y^cos/? M )-2/r + m , (2.1) 

where ^ 

T = -(7i + 72 + 73 + 74) (r 2 = l) (2.2) 

and 

y^=r Yfl r = r-^. (2.3) 

D{p) vanishes at p\ = (0,0,0,0) and p2 = (n/2,n/2,7l/2,7i/2), and can also be seen as a linear 
combination of two physically equivalent naive fermions (one of them translated in momentum 
space). 

The free Karsten-Wilczek Dirac operator is given in momentum space by 

4 3 

D (P) = i L Yti sin/v + i'74 £ (1 - cosp k ), (2.4) 

11=1 k=\ 

'For recent developments, see also jlo|. 

2 They remain rather convenient also for Nf = 2 + 1 and Nf = 2+1 + 1 simulations. The second doublet of minimally 
doubled quarks will contain chirality-breaking terms in order to give different masses to the s and c quarks, however this 
is not so important for these larger masses. 
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and its zeros are instead at p\ = (0,0,0,0) and P2 = (0,0,0, ft). 

The two zeros of these actions, corresponding to the physical flavours, select a special direction 
in euclidean spacetime, identified by the line that connects them. It is easy to see that in the Borici- 
Creutz case the matrix Y selects as a special direction the major hypercube diagonal, while in the 
Karsten-Wilczek case is the temporal direction which becomes the special one. 

As a consequence, hyper-cubic symmetry is broken, and these actions are symmetric only 
under the subgroup of the hyper-cubic group which preserves (up to a sign) the respective special 
direction. This opens the way to mixings of a new kind under renormalization. One of the main 
aims of our work is the investigation of the mixing patterns that appear in radiative corrections. 
We have elucidated the one-loop structure of these theories, and one of our main results is that 
everything is consistent at the one loop level, and the new mixings are very few. 

We also remark that, although the distance between the two zeros is the same {p\ — p\ = ft 2 ), 
these two realizations of minimally doubled fermions are not equivalent. 



3. Counterterms 



Each of the two actions ( |2.lD and ( |2.4| ) does not contain all possible operators which are in- 
variant under the subgroup of the hyper-cubic group preserving its respective special direction. 
Radiative corrections then generate new contributions whose form is not matched by any term in 
the original bare actions. It becomes necessary to introduce counterterms to the bare actions in 
order to obtain a consistent renormalized theory. Enforcing the consistency requirement will allow 
us to uniquely determine the coefficients of these counterterms. 

One must add to the bare actions all possible counterterms allowed by the remnant symmetries. 
Moreover, counterterms are needed also in the pure gauge part of the actions of minimally doubled 
fermions. The reason for this is that, although at the bare level the breaking of hyper-cubic sym- 
metry happens only in the fermionic parts of the actions, in the renormalized theory it propagates 
(via the interactions between quarks and gluons) also to the pure gauge sector. 

We consider the massless case = 0, and note that chiral symmetry strongly restricts the 
number of possible counterterms. It turns out that there is only one possible dimension-four 
fermionic counterterm, which for Borici-Creufz fermions is written in continuum form as 
T/7rX^ D^y. A possible discretization for it has a form similar to the hopping term in the action: 

C4(go) ^^(^( x ) r ^( x )^( x + a ^)~^( x + a A i ) rC/ pi( x )V / '( x ))- (3J) 



There is also one counterterm of dimension three, 



a 



y7(x)r>(x), (3.2) 



3 It is interesting in this respect to observe that an action which contains doublers will in general select some special 
direction, and hence require counterterms. However, the staggered fermion formulation is very clever, because it rear- 
ranges the 16 spin-flavour components of the 4 doublers on the corners of the unit hypercube. Thanks to this, no special 
direction arises, and thus no extra counterterms are needed for the simulation of staggered fermions. In the case of naive 
fermions the 16 doublers are also uniformly distributed in the Brillouin zone, and hence there is no special direction in 
this case too. 
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which is already present in the bare Borici-Creutz action, but with a fixed coefficient —2/ a. The 
appearance of this counterterm means that in the general renormalized action the coefficient of the 
dimension-three operator must be kept general. 

For Karsten-Wilczek fermions we find a similar situation. The only gauge-invariant fermionic 
counterterm of dimension four is 

W74D4Y, (3-3) 

and a suitable discretization of it is 

d 4 {go) ^(^{x)y A U A {x)y{x + a4)-y(x + a4)y A ul(x)y(x)}. (3.4) 
The counterterm of dimension three is for this action 

i ^M W (x)YM^) (3-5) 

(already present in the bare Karsten-Wilczek action, with a fixed coefficient). 

The rules for the counterterm corrections to fermion propagators, needed for our one-loop cal- 
culations, can be easily derived. For external lines, they are given in momentum space respectively 
by 

-;c 4 (go)r£> V; _^M r (3.6) 

v a 

for Borici-Creutz fermions, and by 

-td 4 (go)74P4, — 74 (3-7) 

for Karsten-Wilczek fermions. 

The gluonic counterterms must be of the form trFF, but with nonconventional choices of the 
indices, reflecting the breaking of the hyper-cubic symmetry. It turns out that there is only one 
purely gluonic counterterm, which for the Borici-Creutz action can be written in continuum form 
as 

cp(go) £trF Ap (x)F pT (x). (3.8) 
Apt 

At one loop this counterterm is relevant only for gluon propagators. Denoting the fixed external 
indices at their ends with /I and v, all possible lattice discretizations of this counterterm give in 
momentum space the same Feynman rule: 



-cp(go) 



(Pli + Pv) £/?A - P 2 ~ <V ( £/?A 

A A 



(3-9) 



Contributions of this kind must be taken into account for a correct renormalization of the vacuum 
polarization (see Section ||). 

In the case of Karsten-Wilczek fermions the counterterm which needs to be introduced can be 
written in continuum form as 

d P {g ) £trF pA (x)F pA (x)5 p4 - (3-10) 

pA 
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The Feynman rule for the insertion of this counterterm in external gluon propagators reads 

-dp{go) [pnPv(8nA + 8 V 4)-8 tlv (p 2 8 tl 48 V 4 + pl)] . (3.11) 

In perturbation theory the coefficients of all counterterms are functions of the coupling which 
start at order gf y We will determine (at one loop) the coefficients of all fermionic and gluonic 
counterterms by requiring that the renormalized self-energy and vacuum polarization, respectively, 
assume their standard form (see Sections § and |||). 

Counterterm interaction vertices are generated as well. However, these vertex insertions are 
at least of order gf r and thus they cannot contribute to the one-loop amplitudes that we study here. 
We also want to emphasize that counterterms not only provide additional Feynman rules for the 
calculation of loop amplitudes. They can also modify Ward identities and hence, in particular, 
contribute additional terms to the conserved currents (see Section ^). 

4. Determination of the fermionic counterterms 

Leaving aside for one moment the counterterms, the quark self-energy of a Borici-Creutz 
fermion is given at one loop by 

p 

L(p,m ) = ij}L 1 (p)+mQL 2 (p)+ci(go)-irYpi l + C2(go)-i-, (4.1) 

CL 



where 4 



loga V + 6.80663 + (!-«)(- \oga 2 p 2 +4.7920K)) 



4 log a 2 p 2 - 29.48729 + (!-«)(- log a 2 / + 5.792010 



(4.2) 
(4.3) 



ci (go) = 1.52766- j^Cp, (4.4) 

c 2 (g ) = 29.54170- jjrjpC F , (4.5) 

with Cf = (N 2 — l)/2N c , and a denotes the gauge parameter in a general covariant gauge. The full 
inverse propagator at one loop can be written (without counterterms) as 

E- 1 (/7,m )= (l-Ei) •{#+mo(l-£2 + £l)-;cir£p M -^r}. (4.6) 

We can only cast the renormalized propagator in the standard form 

Z(p,m ) = —^ , (4.7) 

where the wave-function and quark mass renormalization factors are given by 

1-Li) , Z m = l-fz 2 -EiV (4.8) 



For our calculations we have developed programs written in the algebraic computer language FORM 1131 14] 
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provided that we employ the counterterms to cancel the Lorentz non-invariant factors {c\ and C2). 

The term proportional to c\ can be eliminated by using the dimension-four counterterm, 
T/7T y, while the term proportional to c 2 can be eliminated using the dimension-three coun- 

terterm, \ja yry. This amounts to applying the insertions of eqs. (3.6) and (3.7). We thus deter- 
mine in this way that at one loop, for Borici-Creutz fermions, 



c 3 (go) =29.54170. 



16k 2 



C F + O(g 4 ), c 4 (go) = 1.52766- 



167T 2 



C F + 0{gt). (4.9) 



Things work out very similarly for Karsten-Wilczek fermions. In this case the inverse propa- 
gator at one loop (without counterterms) is 



L l (p,mo) = (l -Li) ■ (V+mo (1-E2+L1) -idi y A p A - — 



74 



(4.10) 



where 



Up) 

di(go) 
d% (go) 



16k 2 

„2 



c F 



16k 2 
-0.12554 



log a V + 9.24089 + (1 - a) ( - log a 2 p 2 +4.792010) 
4 loga 2 /? 2 - 24.36875 + (1 - a) ( - loga 2 /? 2 + 5.792010 



16k 2 

„2 



-29.53230- 



167T 2 



C F , 



C F . 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



By using the appropriate counterterms Y74^a¥ an d 1 /a ty 74 ty the renormalized propagator can 
be written in the standard form. Then, at one loop we obtain 



Wo J 



si 



29.53230- C F + 0(jg%), d 4 (g ) = -0.12554- -^C P + 0(gl). (4.15) 



167T 



16?r 2 



One may expect that the above subtraction procedure can be carried out systematically at every 
order of perturbation theory. After the subtractions via the appropriate counterterms are properly 
taken into account, the extra terms appearing in the self-energy disappear. 



5. Determination of the gluonic counterterms 

Leaving aside for one moment the counterterms, the contribution of the fermionic loops to the 
one-loop vacuum polarization of Borici-Creutz fermions comes out from our calculations as 



n$(p) = I PuPv-s^vp 2 



-|Lc 2 ( -~log/>V + 23.6793 



(5-1) 



(/v + Pv) Y,p* ~ p 2 ~ *V ( 



16k 2 



C 2 • 0.9094, 



6 



Minimally doubled fermions 



Stefano Capitani 



where Tr (t a t b ) = C 2 5 ab . For Karsten-Wilczek fermions the corresponding result is 



r C 2 I --\ogp 2 a 2 + 19.99468 



(5.2) 



i67i 2 n 3 

- [puPv (5 M4 + 5 v4 ) - S„ v (/>* + pj) J Jgj C 2 • 12.69766. 

We notice the appearance of non-standard terms, compared with e.g. Wilson fermions. These 
new terms break hyper-cubic symmetry. It is remarkable that they still satisfy the Ward identity 

p"nj5(p) = o. 

At this stage we can employ the gluonic counterterms, which correspond to the insertions in 
the gluon propagator according to eqs. (3.9) and ( p.l 1[ ), to cancel the hyper-cubic-breaking terms 



in the vacuum polarization. The coefficients of these counterterms are hence determined as 

2 2 
c P (go) = -0.9094- -^C 2 + 0{g A Q ), d P (go) = -12.69766- ^C 2 + 0(g 4 Q ). (5.3) 

It is also very important to remark that no power-divergences (1 /a 2 or 1/a) show up in our 
results for the vacuum polarization. 



6. Conserved currents 

We have also calculated the renormalization of the local Dirac bilinears. We have found that 
no mixings occur for the scalar and pseudoscalar densities and the tensor current. For the vector 
and axial currents instead a mixing can be seen, which is a consequence of the breaking of hyper- 
cubic invariance, and their renormalization factors Zy and Za are thus are not equal to one (for their 
numerical values see Section ||). These local currents are indeed not conserved. Using chiral Ward 
identities we have then derived the expressions of the conserved currents, which are protected from 
renormalization. 

As we have previously remarked, the counterterms influence the expressions of the conserved 
currents. It is easy to see that the counterterm of dimension three does not modify the Ward identi- 
ties, and is irrelevant in this regard. On the contrary, the dimension-four counterterm 

£l j^LL (W(x) Yv u n( x ) V(* + <$) + W + ofi) 7v (x) y(jc)) (6. 1) 

^ il V 

generates new terms in the Ward identities and hence contributes to the conserved currents. The 
conserved axial current for Borici-Creutz fermions in the renormalized theory turns out to have the 
expression 

^ W = \ (ft ^ M + a V) + v(* + a #) (7m - * in) 75 uji (*) W, 

+ C4^o) ^) T y 5 U^x)^{x + a^)+W{x + a^)Ty,Ul{x)^{x)y (6.2) 
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For Karsten-Wilczek fermions, application of the chiral Ward identities gives for the conserved 
axial current 

A l (*) = \ \W(x) (Yn ~ m ( 1 - fyO ) 75 Up (x) y(x + apL) 

+y(x + ajl) (Yn + m (1 - V)) ? 5 u l(*) V(*)J ( 6 - 3 ) 
+ M^O) 74 75 ^4 (x) y(x + aA) + y7(* + aA) y 4 75 E/| (x) y (x) 

The conserved vector currents can be obtained by simply dropping the 75 matrices from the above 
expressions. We remark that the vector current is isospin-singlet, representing the conservation of 
fermion number (as also discussed in Jl5|]). The axial current, however, is a non-singlet because 
the doubled fermions have opposite chirality. All these currents have a very simple structure, which 
involves only nearest-neighbour sites. 

We have computed the renormalization of these point-split currents, and verified that is one. 
As all four cases are very similar, we briefly discuss here the conserved vector current for Borici- 
Creutz fermions, for which the sum of the "standard" diagrams (vertex, sails and operator tadpole, 
without the counterterm) gives 



167T 2 



2 



loga 2 p 2 - 6.80664 + (1 - a) ( log a 2 /? 2 - 4.79202 



+ cf(g )r. (6.4) 



The value of the coefficient of the mixing is Cj v (go) = —1.52766 • j^sCf + 0{g^). 

When one adds to this result the wave-function renormalization (that is, £i (p) of the quark 
self-energy), the term proportional to 7^ is exactly cancelled. The mixing term, proportional to T, 
instead remains, because we have not yet taken into account the counterterm. 

The part of the conserved vector current due to the counterterm corresponds to the last line 
of eq. ( |6.2| ). Its 1-loop contribution is quite easy to compute (since C4 is already of order gfy, and 
is given by C4(go)V. We now note that the value of C4 is already known from the self-energy, 
and numerical inspection shows that C4(go) = —c\ v (go) (within the precision of our integration 
routines). Thus, the T mixing term is finally cancelled. We emphasize that only this particular 
value of C4, determined from the self-energy, does exactly this job. 

We have thus obtained that the renormalization constant of these point-split currents is one, 
which confirms that they are conserved currents. Everything turns out to be consistent at the one 
loop level. 



7. Numerical simulations 

If we use the nearest-neighbour forward covariant derivative V^y( x ) = W^( x ) yr(x+dji) — 
y(x)] and the corresponding backward one V* , we can express the (bare) actions in position space 
in a rather compact form. It then becomes apparent that these two realizations of minimally doubled 
fermions bear a close formal resemblance to Wilson fermions: 

<i son = Ut + -"twA, (7-1) 
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D 



KW 



E7m(v m +v*)-^E v * v * ■ 



u=l 



3 

c 



(7.2) 
(7.3) 



All these three formulations contain a dimension-five operator in the bare action, and so we expect 
leading lattice artefacts to be of order a. However, for minimally doubled fermions these effects 



could numerically be small, if the results of [ 16] are to be believed. 

We will not discuss here how to achieve one-loop (or nonperturbative) order a improvement 
for these theories. The classification of all relevant independent operators could turn out to require 
a lengthy analysis. Notice that additional dimension-5 operators will occur not only in the quark 
sector (e.g., "ty/T^yD^Dvi//), but also in the pure gauge part (e.g., Y^nvX F^ V D^F^ V ). Indeed, when 
Lorentz invariance is broken, the statement that only operators with even dimension can appear in 
the pure gauge action is no longer true. 

We would now like to see what can be learned, from the one-loop calculations that we have 
carried out, regarding the numerical simulations of minimally doubled fermions. These simulations 
will have to employ the complete renormalized actions, including the counterterms. 

The renormalized action for Borici-Creufz fermions in position space contains three countert- 
erms and reads 



1 



s bc = « 4 Ei t- E yK*)(yM+c4(j3)r+Y M )tf M (*M*+a£) 



- y(x + a\i ) (y M + c 4 (fi ) r - iy^ ) & J (x) y (x) 



+ \j/(x) ( m + c 3 (j3 ) l — ) Xj/(x) 



E f l-^RetrP MV ] +c P (/3) EtrFjf p (x)Fjt(x)}, 
H<v \ J ,uvp J 



(7.4) 



where F lat is some lattice discretization of the field-strength tensor. We have here redefined the 
coefficient of the dimension-3 counterterm, using C3 (j8) = — 2 + C3(j8) (which does not vanish at 
tree level). 5 

The renormalized action for Karsten-Wilczek fermions also contains three counterterms and 
reads 



KW 



1 



/l=l 



- \ff(x + aji ) (y M (1 + d A (J3 ) 5^) + i/4 (1 - V) ) u l W V W 



We assume that simulations will be carried out at very small values of mo, so that our analysis of the counterterms, 



which assumes chiral symmetry, is essentially still valid. But note also that in our results of eqs. i. 1 



and 4-.1C , obtained 



for general mo, no new dimension-four terms proportional to this mass appear (apart from the standard one, £2)- Thus, 
at one loop we do not need further counterterms in additions to the three which we have found. This strongly suggests 
that our analysis of the counterterms remains valid even when chiral symmetry is broken. 
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C3O0) 


= -2 + 4\ 


8§ + <f< 






=3+4^+4^+ 


C4O0) 


- c?> 


Sl + c?\ 


?$+...; 


^4 (go) 




cp(go) 


= c? 






dp (go) 


- 4V +4 2) 4+ 



+Z 3 I f 1_ ^ Retr vl (l+dp^V)} ( 7 - 5 ) 

(c?3 (y3 ) = 3 + d 3 (p) has a non-zero value at tree level). 

In perturbation theory the coefficients of the counterterms have the expansions 

(7.6) 
(7.7) 
(7.8) 

The same counterterms also appear at the nonperturbative level, and need to be taken into account 
for a consistent simulation of these fermions. Their nonperturbative determination is one the most 
important task for the near future. This can be achieved using suitable renormalization conditions, 
and it remains to be seen which ones will turn out to be more convenient in practice. 

We have previously seen that in perturbation theory the four-dimensional fermionic counter- 
term is necessary for the proper construction of the conserved currents. Its coefficient, as deter- 
mined from the one-loop self-energy, has exactly the right value for which the conserved currents 
remain unrenormalized. This suggests that one possible nonperturbative determination of C4 (and 
<U) can be accomplished by simulating matrix elements of the (unrenormalized) conserved current, 
and imposing (by tuning the coefficient) that the electric charge is one. 

Another effect of radiative corrections is to move the poles of the quark propagator away from 
their tree-level positions. It is the task of the dimension-three counterterm, for the appropriate value 
of the coefficient C3 (or dj), to bring the two poles back to their original locations. These shifts of 
the poles can introduce oscillations in some hadronic correlation functions as a function of time 
separation (similarly to staggered fermions). Then one possible way to determine C3 (^3) is to tune 
it in appropriately chosen correlation functions until these oscillations are removed. 

Such oscillations, familiar from the staggered formulation, come about since the underlying 
fermion field can create several different species, and these species occur in different regions of 
the Brillouin zone. It would be interesting to explore whether or not these oscillations could be 



cancelled by constructing hadronic operators spread over nearby neighbours [17]. 

It is important to remember that because the two species are of opposite chirality, the naive 
7 5 matrix is physically a flavour non-singlet. The naive on-site pseudoscalar field \jfy 5 \ff can create 
only flavour non-singlet pseudoscalar states. To create the flavour-singlet pseudoscalar meson, 
which gets its mass from the anomaly, one needs to combine fields on nearby sites with appropriate 
phases. 

We would like to stress that the breaking of hyper-cubic symmetry does not generate any 
sign problem for the Monte Carlo generation of configurations. The gauge action is real, and the 
eigenvalues of the Dirac operator come in complex conjugate pairs, so that the fermion determinant 
is always non-negative. 

The purely gluonic counterterm for Borici-Creutz fermions introduces in the renormalized 
action operators of the kind E ■ B, E1E2, B2B3 (and similar). In a hyper-cubic invariant theory, 
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instead, only the standard terms E 2 and B 2 are allowed. Fixing the coefficient cp could then be 
done by measuring (E B), (E1E2), ■ ■ ■, and tuning cp in such a way that one (or more) of these 
expectation values is restored to its proper value pertinent to a hyper-cubic invariant theory, i.e. 
zero. These effects could turn out to be rather small, given that only the fermionic part of the tree- 
level action breaks hyper-cubic symmetry. It could also be that other derived quantities are more 
sensitive to this coefficient, and more suitable for its nonperturbative determination. In general one 
can look for Ward identities in which violations of the standard Lorentz invariant form, as functions 
of cp, occur. 

For Karsten-Wilczek fermions the purely gluonic counterterm introduces an asymmetry be- 
tween Wilson loops containing temporal links relative to those involving spatial links only. One 
could then fix dp by computing a Wilson loop lying entirely in two spatial directions, and then 
equating its result to an ordinary Wilson loop which also has links in the time direction. 

In the end only Monte Carlo simulations will reveal the actual amount of symmetry breaking. 
This could turn out to be large or small depending on the observable considered. One important 
such quantity is the mass splitting of the charged pions relative to the neutral pion. Indeed, since 
there is only a£/(l)®I/(l) chiral symmetry, the 7T° is massless, as the unique Goldstone boson 
(for mo — > 0), but n + and k~ are massive. 

Furthermore, the magnitude of these symmetry-breaking effects could turn out to be substan- 
tially different for Borici-Creutz compared to Karsten-Wilczek fermions. Thus, one of these two 
actions could in this way be raised to become the preferred one for numerical simulations. 

8. A unifying notation for the two fermion discretizations 

By introducing a particular notation, some similarities between the two realizations of min- 
imally doubled fermions can be revealed. This applies to the form of the action, operators and 
counterterms. For this purpose one can introduce a 4-component object A^, defined as 

f 8n4 Karsten-Wilczek f y 4 Karsten-Wilczek 

^ i Borici-Creutz ^ F Borici-Creutz 

In both cases this object points from the zero of the action at the center of the Brillouin zone to the 
other zero (describing the second fermion, of opposite chirality). 

At first we show that by means of this object one can cast both actions into similar (although 
non-equivalent) forms. Their free Dirac operators, as we have already seen in Section [7J, contain 
the same naive fermion piece but a different dimension-five operator. The latter can be rewritten in 
this new notation as 

CI „ y Z \ / 

J>2k*) - --lAVsin 2 ^(l-25, v ). (8.3) 

The factors (1 — 8^ v ) and (1 —28^) cannot be transformed into each other, and this illustrates that 
the two actions are inequivalent and must be distinguished (as we remarked in Section ||). 
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Although the quark propagator cannot be cast into a uniform expression using this notation, 
this turns out to be possible for operators (e.g. local currents and counterterms), as well as some 
other results such as the expression for vacuum polarization. For example, the various counterterms 
that we have previously discussed can be easily cast in a completely unified way for the two actions. 
If we rewrite the three counterterms making use of the object A fl , the counterterms of dimension 
three appear as 

iW(x){A-r)y(x), (8.4) 
the fermionic ones of dimension four become 

W(x)(A-y)(A-D) ¥ (x), (8.5) 

and the gluonic ones are 

£ A^ p F pv A v . (8.6) 

jU,V,p 

Here (and in the following) objects written in this unified notation may differ by simple numerical 
coefficients from the corresponding quantities which we have previously used in the conventional 
notation. 

Let us now consider the results of the one-loop calculation that we have presented in the 
previous Sections. One can rewrite the full self-energy (without counterterms) for both actions as 

E = ip 1 Ej + m E 2 + ici (A • y) (A • p) + c 2 (A • y) , (8.7) 

with di being given by either c, or d\. Also the fermionic bilinears can be expressed in a unified form. 
Using the abbreviations b = and L = loga 2 p 2 , the results for the one-loop vertex diagram for 
the local scalar, vector and tensor bilinears are 

- 4L + 24.36875 + (1 - a) (L- 5.792010) ) Karsten-Wilczek 
(y> — j v v ' ) (8 8) 

1 (- 4L + 29.48729 + (1- a) (L- 5.792010)) Borici-Creutz 
y,(-L + 10.44610+ (l-a)(L-4.792010)) - 2.88914 -A M (A -y) Karsten-Wilczek 



L + 9.54612 + (1- a) (L- 4.792010)) - 0.20074 -A M (A -y) Bori^-Creutz 

a uv f4.17551 + (l - a) (1-3.792010)") Karsten-Wilczek 

> < • (8-10) 

a MV ( 2. 16548 + (1 - a) (L- 3.792010) I Borifi-Creutz 



For the conserved vector current, the sum of the standard proper diagrams (vertex, sails and opera- 
tor tadpole) reads for the two actions 



b I 



- L - 9.24089 + (1 - a) (L - 4.792010) + 0. 12554 ■ A M (A ■ y) Karsten-Wilczek 

- L - 6.80663 + (1 - a) (L - 4.792010) ) - 3.05532 • A^ (A • y) Borici-Creutz 

(8.11) 

Perhaps one of the most striking examples of the convenience of this notation can be observed 
in the case of the vacuum polarization. The contribution of fermion loops to this quantity contains 
structures which break hyper-cubic symmetry. It can be written as 

11$ (p) = I 3 (PuPv - P 2 8^ v ) + d g ((A • p) (A^ Pv + A vPll ) - (A^Avp 2 + 5 MV (A-/?) 2 )), (8. 12) 
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with the numerical results (as we have seen in Section 



2\ _t) ~\ L + 19.99468 Karsten-Wilczek 



SO 



3 



L + 23.6793 Borici-Creutz 



, , r -12.69766 Karsten-Wilczek 

= " -3.6376 B„n si -Creu«z ' (8 '' 4) 



with = (Wilson fermions have £3 = ft(—|L + 4.337002) and cfg = 0). Thus, a single for- 
mula can describe the structures which arise in the calculation of the vacuum polarization for both 
actions. 

With this notation we have thus shown that operator structures and results for Borici-Creutz 
and Karsten-Wilczek fermions, although distinct, share many common traits. As can be inspected 
in the above expressions, another remarkable feature appears to be that, after is introduced, the 
summed indices occur in pairs (like in the continuum), and also the free indices match exactly on 
both sides of equations. We do not know if this will always happen, also if one computes more 
complicated quantities. 

Even without using the A notation, we also discovered that the hyper-cubic-breaking terms 



of the vacuum polarization in eq. (8.12) can be put for both actions in the same algebraic form, 
namely 

p 2 {r^n{Yv,n + s^,r}{^r}-^j}^j}{Y v ,r} + { Yv j}{Y^r}), (8.15) 

where in the case of Karsten-Wilczek fermions T must be replaced by 74/2. This substitution is 
suggested by comparison of the standard relation F = | L^i(7m + 7u) of Borici-Creutz fermions 
with the formula 74 = 5 L^i(7m + f° r Karsten-Wilczek fermions, expressing the symmetries 
of the action (as can be seen from the expression of the propagator, when one expands it around 
the second zero). Whether there is any deeper significance to this structural "equivalence" of the 
hyper-cubic-breaking structures in the vacuum polarizations remains an open question. 



9. Conclusions 



Borici-Creutz and Karsten-Wilczek fermions are described by a fully consistent renormalized 
quantum field theory. Three counterterms need to be added to the bare actions, and all their co- 
efficients can be calculated either in perturbation theory (as we have shown), or nonperturbatively 
from Monte Carlo simulations (a task for the future, for which we have suggested some strategies). 
After these subtractions are consistently taken into account, the power divergence in the self-energy 
is eliminated, and no other power divergences occur for all quantities that we calculated. 

We have argued that under reasonable assumptions and following the nonperturbative deter- 
mination of these counterterms, no special features of these two realizations of minimally doubled 
fermions should hinder their successful Monte Carlo simulation. 

Conserved vector and axial currents can be derived, and they have simple expressions which 
involve only nearest-neighbours sites. We have then here one of the very few cases where one can 
define a simple conserved axial current (also ultralocal). 
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Finally, we would like to observe that this work is also an example of the usefulness of pertur- 
bation theory in helping to unfold theoretical aspects of (new) lattice formulations. 
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